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We have obtained strong experimental evidence for the full determination of the superconducting 
gap structure in all three bands of the spin-triplet superconductor Sr2Ru04 for the first time. We 
have extended the measurements of the field-orientation dependent specific heat to include conical 
field rotations consisting of in-plane azimuthal angle (ji-sweeps at various polar angles 6 performed 
down to 0.1 K. Clear 4-fold oscillations of the specific heat and a rapid suppression of it by changing 
9 are explained only by a compensation from two types of bands with anti-phase gap anisotropics 
with each other. The results indicate that the active band, responsible for the superconducting 
instability, is the 7-band with the lines of gap minima along the [100] directions, and the passive 
band is the a- and /3-bands with the lines of gap minima or zeros along the [110] directions in their 
induced superconducting gaps. We also demonstrated the scaling of the specific heat for the field 
in the c-direction, which supports the line-node-like gap structures running along the kz direction. 



I. INTRODUCTION 



The unconventional superconductivity, realized in 
high-Tc cuprates, organics, heavy-fermion intermetallic 
compounds, ruthenate, etc., has become one of the most 
actively studied topics of condensed matter physics to- 
day. Since the discovery of its superconductivity,^ the 
layered ruthenate Sr2Ru04 has generated extensive re- 
search activities. What make Sr2Ru04 special among 
many unconventional superconductors are the availabili- 
ties of detailed knowledge of its Fermi-liquid properties, 
as well as of ultra-high-purity single crystals with the 
residual resistivity less than 100 nficm.'* In particular, 
its Fermi surface is relatively simple and consists of three 
nearly-cylindrical sheets of a-, /?-, and 7-bands with spe- 
cific d-electron orbital characters. Fermi-surface param- 
eters of all three sheets, reflecting strong electron corre- 
lations, are experimentally quantified in detail.^'^ There- 
fore, there is a prospect of the underlying physics of its 
unconventional superconductivity be understood in un- 
precedented detail and depth. Obviously, the determi- 
nation of the symmetry and anisotropic structure of the 
order parameter, the superconducting gap, is a crucial 
step toward identification of its pairing mechanism. 

The superconductivity of Sr2Ru04 has pronounced un- 
conventional features such as: the invariance of the spin 
susceptibility across its superconducting (SC) transition 
temperature Tc,"^'^ appearance of spontaneous internal 
field, ^ evidence for two-component order parameter from 
field distribution in the vortex lattice, absence of a 
Hebel-Slichter peak,^^ and the strong dependence of Tc 
on nonmagnetic impurities. "'^^'^'^ These features are co- 
herently understood in terms of spin-triplet p-wave su- 
perconductivity with the vector order parameter being 
d{k) = zAo{kx H- iky),^'^ in contrast to the spin-singlet 
rf-wave pairing known in high Tc cuprates. This order 
parameter represents the paired spin state with S ~ 1 
and Sz = 0, and the orbital wave function with = +1, 



called the chiral p-wave state. The above simple vector 
order parameter leads to the gap A{k) = Ao{kl-\- ky)^^'^ , 
which is isotropic because of the quasi-two dimensional- 
ity of the Fermi surface. Exponential temperature de- 
pendence of quasiparticle (QP) excitations is expected in 
this case. 

However, a number of experimental results with very 
high quality crystals revealed the power-law temperature 
dependence of QP excitations. In particular, the specific 
heat^^ C/T cx T , the thermal conductivity^^ k/T cx T 
and NQR relaxation rate^^ l/Ti c>c T^, in addition to 
the ultrasonic attenuation rate^^ a and variation of the 
penetration depth-"^^ AX, all indicated the presence of the 
lines of nodes, of zeros, or of gap minima in the supercon- 
ducting gap. Nevertheless, magnetothermal conductivity 
measurements for determination of the SC gap structure, 
with the applied field rotated within the RUO2 plane 
down to 0.35 K, have indicated little anisotropy within 
the planes. In spite of a number of detailed studies 
of the SC gap structure, the wave function of Sr2Ru04 
are thus still controversial. 

The main outstanding issues on the unconventional 
superconductivity of Sr2Ru04 at present may be sum- 
marized as follows^: (1) full characterization of the SC 
gap structure (i.e. orbital state) in all bands, (2) di- 
rect observation of odd parity of the orbital state by a 
phase-sensitive method (3) clarification of the details of 
the spin state in the spin-triplet pairing state especially 
under magnetic fields, (4) confirmation of the "chiral- 
ity" due to time reversal symmetry breaking in the SC 
state and searches for novel phenomena associated with 
chirality, and (5) understanding of the multi-phase su- 
perconducting phases and the critical field limiting. ^^"^^ 
Furthermore, the assignment of the dominant band in 
producing the superconductivity by taking the orbital 
dependent superconductivity (ODS) into account^'^, as 
well as the examination of the role of the incommensu- 
rate AF fluctuation to pairing in the triplet channel, is 
crucial for the determination of the pairing mechanism. 
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We report in this paper a new development to resolve 
the issue (1). Measurements of the spectra of QPs in- 
duced by magnetic fields serve as a powerful approach 
for this purpose. The field-orientation dependent spe- 
cific heat is a direct measure of the QP density of states 
(DOS) and thus a powerful probe of the SC gap struc- 
ture.^*'"'^" Indeed, in our previous report, we revealed the 
emergence of in-plane anisotropy below about 0.3 K in 
the SC gap of Sr2Ru04 from the field-orientation depen- 
dence of the specific heat.^^ We identified the multi-gap 
superconductivity with the active band 7, which has a 
modulated SC gap with a minimum along the [100] di- 
rection. 

In the present study, we extended the measurements 

of the field-orientation dependence of the specific heat to 
the field direction covering all polar and azimuthal an- 
gles. We also found the scaling of electronic specific heat 
in magnetic fields perpendicular to the RUO2 plane. As 
a result, we have determined the superconducting gap 
structure of the spin-triplet superconductor Sr2Ru04 in 
all bands for the first time. We will examine and se- 
lect the pertinent theoretical models compatible with the 
present observation and discuss the effective pairing in- 
teraction as well as the roles of the incommensurate AF 
fluctuation in a- and /3-bands. 



II. EXPERIMENT 

Single crystals of Sr2Ru04 were grown by a floating- 
zone method in an infrared image furnace.^ After speciflc- 
heat measurements on two crystals, the sample with 
Tc = 1.48 K, close to the estimated value for impurity 
and defect free specimen (Tco = 1.50 K),^^'^^ was cho- 
sen for detailed study. This crystal was cut and cleaved 
from the single crystalline rod, to a size of 2.8 x 4.8 mm^ 
in the a6-plane and 0.50 mm along the c-axis. An x-ray 
rocking curve of the sample shows the characteristics of a 
single crystal of high quality; the diffraction peak width 
[full width at half maximum (FWHM)] was comparable 
to that of a Si crystal (with FWHM of 0.06°) in our 
diffractometer. The directions of the tetragonal crystal- 
lographic axes of the sample were determined by x-ray 
Laue pictures. The side of the crystal was intentionally 
misaligned from the [110] axis by 16°. 

The field-orientation dependence of the specific heat 
was measured by a relaxation method with a dilution 
refrigerator. The details of the measurement system, 
consisting of two orthogonally-arranged SC magnets and 
a mechanical rotating stage for the dewar, are described 
elsewhere. "^^ In the present study, we mainly operate the 
system with constant polar angle 9 and varing azimuthal 
angle cf) with a misalignment AO no greater than 0.01°. 
It is important to maintain a constant 9 during this 
conical variation of the field direction, because the 
large ifc2 anisotropy in 6. The in-plane and inter-plane 
field alignment relative to the crystallographic axis was 
carried out, based on the known in-plane anisotropy 



of Hc2 at low temperatures. ^^"^^ For referring to the 
direction of the applied magnetic field with respect to 
the crystallographic axes, we shall introduce the polar 
angle 9, for which = 0° corresponds to the [001] 
direction, and the azimuthal angle (j), for which = 0° 
with 9 = 90° corresponds to the [100] direction. It is to 
be noted that the tetragonal symmetry of the crystal 
structure is known to conserve down to temperatures as 
low as 110 mK.33 



III. RESULTS AND DISCUSSION 
A. Scaling of the specific heat for H || [001] 

The objective of this section is to present the QP 

DOS in the vortex state of Sr2Ru04, when a magnetic 
field is applied perpendicular to the conducting plane, 
H II [001]. Prom such information, we try to determine 
the dimension of nodes or of gap minima: lines or points. 
Figure 1(a) shows the temperature dependence of the 
electronic specific heat divided by temperature Ce/T at 
several fields. The electronic specific heat Ce in mag- 
netic fields was obtained after subtraction of the phonon 
contribution with a Debye temperature of 410 K. The 
temperature dependence of Ce/T at low temperatures is 
proportional to T down to 50 mK, and this tempera- 
ture dependence suggests lines of nodes or of gap min- 
ima in the SC gap.^^ Note that Cc/T at zero field lin- 
early extrapolates to the residual electronic specific-heat 
coefficient 70 from about 0.5 K with Cc/T = 70 + oT: 
70 = 3.8 mJ/K^mol and a = 48.8 mJ/K'''mol. The re- 
sult of 70/7N = 0.1 guarantees that this sample is clean 
enough, where 7n is the electronic specific heat coefficient 
in the normal state: 7n = 37.8 mJ/K^mol. 

To identify the dominant source of the DOS at low 
temperatures, we estimate the impurity scattering rate 
divided by the SC gap amplitude ftr/Z\o = 0.0098 
for Tc = 1.48 K from the Abrikosov-Gorkov equation 
ln(Tco/Tc) = ^(1/2 + hr/2nkBT,) - *(l/2). * is a 
digamma function, Tco = 1.50 K is a maximum Tc 
for the disorder- free material, and we approximate 
2Z\o/A:bTco — 3.53 (BCS). This result is in good agree- 
ment with the residual DOS estimated from the unitarity 
limit: 70/7N — ^/h^ / Aq.^^'^^ The condition of the su- 
perclean regime <C < can be satisfied above 
Hci, where Hc2o ^ ^fc2(T = o'k) = 71 mT for H \\ c 
is estimated from ffc2(T) determined from the specific 
heat in Fig. 1(a) and (b). Thus we conclude that the 
effect of impurity scattering is negligible. For < , 
QPs are dominated by those induced by magnetic fields; 
in the mixed state, the QP energy spectrum Ek is af- 
fected and becomes E'j^ = Ep, — Slo by the Doppler shift 
Soj = hk ■ Vs, where Vg is the superfluid velocity around 
the vortices and hk is the QP momentum. In the case 
of Suj > A{k), this energy shift gives rise to a finite DOS 
at the Fermi level. This leads to the field dependence of 
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FIG. 1: (a) Temperature dependence of the electronic specific 
heat divided by temperature Cc/T at several fields for H || 
[001]. The line at zero field is a fit of the low temperature 
Cc/T with 7o + aT. (b) Field dependence of Cc/T at several 
temperatures for H \\ [001]. The line is a fit with 70 + AsfH 
for the data at 60 mK . 



Cc/T cx DOS being proportional to at low tempera- 
ture for the SC gap with lines of nodes or of gap minima 
(i.e. for the QP DOS of cx ^, where E is the en- 
ergy of a QP)^^: (X at zero temperature, where 

Eh — Ao^H/Hc2o is the energy scale associated with 
the Doppler shift Suj.^^ Np is DOS at the Fermi level in 
the normal state. 

Figure 1(b) shows the field dependence of Cc/T at sev- 
eral temperatures for H || [GDI], measured after zero-field 
cooling. At low temperatures, the specific heat behaves 
like Ce/T (X ^/H. The deviation from the Vh law for 
H < 14 mT is most likely due to entrance into the Meiss- 
ner state (iJci(O) ~ 7 mT). Well above Hd, field induced 
QPs become dominant and C'e/T has little temperature 
dependence below 0.1 K as shown in Fig. 1(b). If the 
data at 60 mK is chosen as the limit of T/Tc 0, the 
field dependence of Ce/T is described as 70 + A^/H with 
A = 4.3 mJ/K^T^/^mol. This field dependence also in- 
dicates nodes or of gap minima in the SC gap. 

Recently, it has been shown that the field induced 
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FIG. 2: (a) Experimental data of Co/T for H \\ [001], 
as a function of field strength and temperature. A con- 
tour plot is shown on the bottom H — T plane, with the 
same color scale as the 3D plot. (b) Parameter x = 
a{T /Tc) / ^ H/Hc20 dependence of the normalized QP con- 
tribution Cqp{T,H)/Cqp{T,0) for H \\ [001] with the same 
color as Fig. 1(a) and (b), where Cqp{T, H) is defined as 
Ca{T,H) — -yoT. The inset shows H-T phase diagram for 
H II [001]. The scaling of Cqp{T, H)/Cqp{T,0) is carried 
out in the shaded region. 



QP DOS should yield general scaling relations of vari- 
ous thermodynamic and transport properties. '^'''^^ When 
the SC gap has lines of nodes or of gap minima in the 
field direction (i.e. lines of nodes or of gap minima || c- 
axis and H \\ c-axis), the QP DOS is a simple func- 
tion of E/y/H. In this case, the scaling function has 
been explicitly calculated for the specific heat. The scal- 
ing relations are functions of a single reduced parameter 

2; — _2_ , T _ „ T 

v-p 



T _ 



■^^ff^, where 0o is the flux quan- 
tum, vp is the Fermi velocity in the conducting plane 
and a is a constant of the order of unity. '^^ Experimen- 
tally, this scaling of the specific heat has been demon- 
strated in high Tc cuprate superconductors.'^^ To focus 
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on the QP excitation only by the temperature and mag- 
netic field, we define CQp(r, H) = C^{T, H)-joT. Then 
Cqp(T, H)/Cqp{T, 0) would be expressed by the scaling 
function F{x) for the SC gap with lines of nodes or of 
gap minima, where F{x) is well approximated with the 
empirical interpolating function i^int(a;) 
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nodes, of zeros, or of gap minima. In addition, the ob- 
servation of Ce/T oc ^/H up to near Hc2 is also consis- 
tent with the chiral p-wave state. *° Note that although it 
might be possible to explain these features by the single 
band model with lines of nodes or of gap minima, it is 
not so simple because the features of multi-band super- 
conductivity emerges only by tilting the field direction to 
near a6-plane, as discussed in next section. 
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F{x) and i^int(a^) has two regimes, depending on whether 
the thermal energy ^^f" — ^ is small {F{x), Fiy^t{x) 
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J for x < 1, Cqp = ATVH for T ^ 0) or large 
(F(x),Fint(a:) ^ 1 for a; » 1, Cqp = aT^ in zero field) 
compared to the energy scale of the Doppler shift ^ ~ 

The asymptotic form of Cqp (T, H) is 



Eh _ _ 
^0 ~ V 

given by 



H 



B. Field-orientation dependence 
of the specific heat 

In the previous section, we establish the scaling of 
Cc/T with the single parameter x oc T/^/H for H \\ 
[001], which indicates the existence of lines of nodes or 
of gap minima in the SC gap. In this section, we present 
the field-orientation dependent specific heat, from which 
we try to determine the direction of lines of nodes or of 
gap minima. Let us first present the results for 6 ~ 90°^^ 
with additional analysis. We will then present the new 
results for 6* 7^ 90°. 

Figure 3(a) shows the field dependence of Ce/T for 
H II a5-plane at T = 0.12 K. Cc/T increases sharply 
up to the characteristic field /io-f^c2 — 0-15 T and then 
increases more slowly for higher field. This behavior is 
more pronounced at 60 mK,"^^ and naturally explained 
by the presence of two kinds of the gap.^^ In view of 



CQp(r,iJ) = CQp(r,0)^^(x) ^ CQp(T,0)Fi„t(x) (4) 
'^27C^)'^\/ TT^ (a; < 1, ^ < ^) 

, (5) 

aT^ (a; > 1^ ^ < 1) 

where CQp(r, 0) approaches aT^ for ^ < i in Sr2Ru04. 

Figure 2(a) is a 3-dimensional (3D) representation of 
Ce/T for the [001] field direction, as a function of field 
and temperature. The figure is constructed from data in 
Fig. 1(a) and (b). Provided that the SC gap structure 
of Sr2Ru04has lines of nodes or of gap minima, but not 
points of those, Cqp (T, H) / Cqp (T, 0) , deduced from this 
full map of the specific heat, would well be described by 
the scaling function F{x) with the reduced parameter 
X = aT/Tc^Hc2a/H in the range of < T < Tc and 

Hel{T)<H <He2{T). 

Figure 2(b) shows the x dependence of the normalized 
QP contribution Cqp(T, ir)/CQp(T, 0) of the electronic 
specific heat for H \\ [001]. Indeed the data obtained 
by i?-sweep at constant T and T-sweep at constant H 
collapse onto a single curve in a wide range in the H- 
T phase diagram, shown in the inset of Fig. 2(b) as the 
shaded region, and the curve is well described by F{x) 
even in the crossover region a; « 1. The constant a is 
determined from a|f^ = A = 4.3 mJ/K^T^/^mol in 

the limit of a; <C 1, where C(3) is 1.202..., and results in 
a = 0.8, indeed an order of unity. Such demonstration 
of the scaling strongly supports the existence of lines of 




(a) • H // [1 1 0] 

■ OH// [100] 



^5 45 90 1350 

^(deg.) 




0.5 1.0 1.5 

/^H{T) 



FIG. 3; (a) Field dependence of Ce/T at T = 0.12 K for 
H II [100] (open circle) and if || [110] (closed circle), (b) 
The in-plane field-orientation dependence of C^/T at several 
fields. The solid lines are fits with /4(0) given in the text, 
(c) Field dependence of the normalized 4-fold oscillation am- 
plitude C4/CN. The points are evaluated from the fitting to 
the oscillatory data in Fig. 3(b), while the lines from the dif- 
ference in Co between H || [110] and H || [100] in Fig. 3(a). 
Two methods yield consistent results. 
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the QP excitations in ODS,^^ both active and passive 
bands are effective below ^-t ^ow temperature region, 
but only the active band is important above and/or 
high temperatures since most of the QPs have already 
been excited for the passive bands. On the basis of the 
different orbital characters of the three Fermi surfaces 
(a, /3, and 7),^'^ the amplitudes of the gap functions Aai3 
and are expected to be significantly different and the 
normalized DOS of those bands are = 0.43 and ^ = 
0.57. The plateau on the field dependence of Cc/T at low 
temperature divides 7n (oc A^f) reasonably into 0.43 x 7n 
recovered at low field and 0.57 x 7n recovered at high 
field. Consequently we conclude that the active band 
which retains strong superconductivity up to high field 
is the 7-band mainly derived from the in-plane dxy orbital 
of Ru M electrons. 

As mentioned in the previous section, the QP energy 
spectrum is affected by the Doppler shift in the mixed 
state for the superfluid velocity around the vortices, and 
consequently this energy shift gives rise to a finite DOS 
at the Fermi level where the local energy gap becomes 
smaller than the amplitude of the Doppler shift. '^^ Since 
the superfluid velocity Vs is perpendicular to the mag- 
netic field H, the Doppler shift Suj becomes zero for 
k II H. Thus the generation of QPs by magnetic fields 
at nodes or at gap minima is suppressed for || node or 
gap minima directions, and Cg/T takes minima for these 
directions. 

Figure 3(b) shows the field-orientation dependence of 
Ce/T for H II a&-plane at T = 0.12 K. The absence of a 
2-fold oscillatory component in the raw data at ^ = 90° 
guarantees that the in-plane field alignment is accurate 
during the azimuthal-angle </> rotation. For the field range 
0.15 T < i^oH < 1.2 T, where the QPs in the active 
band 7 are the dominant source of in-plane anisotropy of 
Cc{(t>), a non- sinusoidal 4-fold angular variation approxi- 
mated as Co(0) = Ca + Cifi{(f)) with f^{(f)) = 2|sin2(/)| - 1 
is observed. Since we observe (i) no angular variation 
above Hc2 in the normal state and (ii) a 4-fold angular 
variation with the phase inverse to that due to the in- 
plane Hc2 anisotropy below ifc2 in the SC state, we con- 
clude that the 4-fold oscillations at low temperature orig- 
inate from the anisotropy in the SC gap.^'^ In addition, 
cusp- like features at the minima {(j) = f n, n: integer) are 
attributable to the strong reduction of QP excitations for 
the nodal direction parallel to H since the Doppler shift 
becomes 5uj — Q for a purely two-dimensional gap struc- 
ture.^^ On the other hand, strong kz dependence of the 
gap function would enhance the QP excitations even for 
the nodal direction parallel to H, so that the cusp-like 
features would have been strongly suppressed. There- 
fore, in the active band 7, the gap structure with vertical 
line nodes or gap minima along the [100] directions are 
promising. In addition, from the consistency with ob- 
served value of Z\Ce/7N21, which mainly originates from 
the active band 7, we can deduce the order parameter 
d{k) = zZ\o(sinafca; -|- isinafcy) with gap minima for the 
active band.^^ 



Figure 3(c) shows the field dependence of C4/CN for 
H II a6-plane at T — 0.12 K. C4/CN is the normal- 
ized angular variation amplitude, where Cn = ^nT is 
the electronic specific heat in the normal state. At low 
fields {jjLoH < 0.15 T) and low temperatures (T < 0.3 K) 
where QPs on both the active and passive bands are im- 
portant, the oscillation anisotropy rapidly decreases. We 
define the threshold field ffth for the strong reduction 
of C4/CN in Fig. 3(c). While the previous experimen- 
tal study^^ has resolved the directions of gap minima in 
the active band 7, there still remains two types of pos- 
sibilities predicted for the passive bands a and /?: (A) 
horizontal line nodes''^"'*^ and (B) vertical lines of nodes 
or of gap minima^^"^^. 

In view of the gap structure in the active band, this 
steep reduction of C4 may be explained with the gap min- 
ima z4niin of d{k) = £Z\o(sinafca; -I- \sinaky): at low field 
(Doppler shift 5uj < /imin) the 4-fold oscillations cannot 
occur, whereas above the threshold field Hth {Suj > Zimin) 
the 4-fold oscillations will be observed. In this case, the 
gap structure of both (A) and (B) are still possible in 
passive bands, since these features may be explained with 
only the active band except for the power-law dependence 
of Ce/T (X T. Note that according to Fig. 3 it is assumed 
in this scenario that H^h characterizing the gap minima 
of the active band happens to be nearly the same as 
characterizing the passive band. 

In view of ODS, the existence of vertical lines of nodes 
or of gap minima in the passive bands may give rise to 
4-fold oscillations in Cc(0) which are out of phase with 
the 4-fold oscillations originating from the active band. 
In this case, the 4-fold oscillations would be suppressed 
below the characteristic field -ff^2j where both active and 
passive bands contribute to the Cc{(p) oscillation. This 
idea requires strong in-plane anisotropy of the gap in pas- 
sive bands and lead to the theoretical models (B). 

We have extended the specific heat measurements to 
cover all 9 and ^, in order to determine the gap struc- 
ture in passive bands. Now, let us focus on the rela- 
tion between the characteristic field Jf^2 of ODS and 
the threshold field i?th of steep C4 reduction. First, 
we trace the change of H[.2 by tilting the polar angle 
9. Next, we investigate the azimuthal angle depen- 
dence of Co at several polar angle 0, fixing the tem- 
perature and the reduced magnetic field H/Hc2{9)- In 
particular, we kept the reduced field such that iJth is 
always exceeded: ^ > ^f^^ ~ ^ i.e. -n-^ > 

W) ^ ti^ = = 0-1' wh ere averaged 

can be estimated from ^ ^ ^ = ^J H / HcioiO) and 
V Hih{9) / Hc2o{9) ^ ^ ^ = constant in the case 
that Hth is determined only by gap minima. ^^"^^ In the 
case that C4 reduction is caused solely by gap minima in 
the active band, the polar angle 9 dependence of C4 is 
expected to be similar to that of single band models with 
vertical lines of nodes under the condition of 5co > /:imin- 
In that case C4 (9) is not expected to be strongly affected 
by the evolution of the characteristic field H'^2 of ODS. 
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FIG. 4: (a) Field dependence of Ce/T at several polar angle 
from the [001] to [100] directions, fixing the temperature 
T = 0.10 K and the azimuthal angle = 0°. (b) Cc/T vs 
^H/Hc2{9)at several polar angle 6 from the [001] to [100] 
directions. 



On the other hand, in the case that C4 reduction is caused 
mainly by compensation of anti-phase gap anisotropics 
between active and passive bands, C4{0) is expected to 
be deeply related to the evolution of Especially, 
C4 will rapidly be suppressed for H < H[.2, for which 
anisotropic QP excitations can occur in both bands. 

First, to trace the variation of with 0, we inves- 
tigate the field dependence of Ce/T by tilting the po- 
lar angle 9. Figure 4(a) shows the field dependence 
of Cc/T at several polar angles 9 from [100] to [001] 
directions, fixing the temperature T = 0.10 K. For 
H II [100], Cc/T increases sharply up to the characteris- 
tic field = 0.15 T and then increases more slowly 
for higher field, which is naturally explained by a pres- 
ence of two kinds of the gap. By tilting the polar angle 
9 from a5-plane, these features gradually change in the 
range of 60° < 9 < 90°, and finally, the field depen- 
dence of Cc/T approach Cc/T oc \f}i in the range of 
0° < < 60°. To extract the polar angle dependence 
of i?c2(^)j we plot in Fig. 4(b) the field dependence of 



Cc/T against yH/H^iW) several polar angle 9 from 
[100] to [001] directions, fixing the temperature T = 0.10 
K. We directly find the variation of iJ^jC^); niarked with 



0.1 .: 



O 
CT) 



1.0 



0.5 



(a) • H, 



c2 
c2 



r= 0.10 K 



▼ measurement 
points of Cg((^ ) 



H, 



c2 



7 h; 



c2 



[001] 



HIH^2 = 0-4 



[100] 



(b) 



3D - line nodes 



r=0.12K 
H/H^2 = 0.4 



/2D - line nodes / 

— g "' i — ^ 




30 60 

^(deg.) 



FIG. 5: (a) Polar angle 6 dependence of Hc2 and at 0.10 
K. The red and blue lines represent fits of Hc2 and Hc2 with 
the Ginzburg-Landau anisotropic effective mass approxima- 
tion. The black line indicates H/Hc2 = 0.4, for which con- 
ical azimuthal angle 4> dependence of Cc/T is investigated, 
(b) Polar angle 6 dependence of the normalized 4-fold oscil- 
lation amplitude dlfi) /C4,{'dQ°) and C4,{6)/Cn in the specific 
heat. The dotted line and dashed line are calculation results 
of C4(6i)/C4(90°). 



arrows in Fig. 4(b), toward the Hc2{9) by tilting the po- 
lar angle 9 from afo-plane in the range of 90° to 60°: 
H'c2i^) ~^ Hc2{9). It is not easy in this way to precisely 
define Hl,2{(^)] we will confirm the validity of this defini- 
tion as described in the next paragraph. In the range of 
Q° < 9 < 60°, we are able to find good scaling such that 
Cc/T shows universal H / Hc2{9) dependence. In this 
region, we judge H'^2i^) to be corresponding to Hc2{9)- 
H'c2{9) = Hc2{9). 

Figure 5(a) shows the polar angle 9 dependence of Hc2 
and H'^2^ determined from the field dependence of the 
specific heat at 0.10 K. For orientations not very close to 
H II a6-plane, the angle dependence of Hc2 is well fitted 
with the Ginzburg-Landau anisotropic effective mass ap- 
proximation."*^ Also shown in Fig. 5(a) is the curve based 
on the G-L effective mass model, 

Tic2\\c 



Hc2{e) 



Vcos2 9 + Ta"^ sin^ 9 



(6) 



Here Ta is the square root of the ratio of the effective 
mass for in-plane motion to that for inter-plane motion 
(i.e. Ta = Hc2\\ab/Hc2\\c)-^^ Here we have taken only 
Hc2\\ab to be the adjustable parameter for the fitting 
with Hc2\\c — 0.071 T given. In reality, Hc2 becomes 
nearly independent of 9 and strongly undershoots the 
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FIG. 6: The azimuthal angle (p dependence of Ce/T at several 
polar angles by rotating H conically, fixing the temperature 
T = 0.12 K and the reduced magnetic field H/Hc2 = 0.4. The 
solid lines are fits with f4{4>) given in the text. 



value obtained from extrapolation of the fitting in the 
88° < 6 < 90° range, due to the hmiting of Hc2 and the 
occurrence of the second superconducting transition in 
the H II Ru02 plane. The effective upper critical field 
of the passive band can also be fitted without any fit 
ting parameter, using F'^ = H'^2\\ab/K2\\c 



H\u . = 0.15 

c2\\ab 



T, and i7^2||c = Hc2\\c = 0.071 T. The fitting curve in 
Fig. 5(a) indeed confirms that -ffc2(^) defined in Fig. 4(b) 
is in accord with the G-L effective mass model. 

The in-plane and intcr-plane coherence lengths, 
and are given in the anisotropic effective mass ap- 
proximation: Hc2\\c = Hc2\\ab = 2ntat^^ " ^^^^^ 
specific heat for || c does not indicate a presence of two 
kinds of the gap, it is reasonable to assume that ifc2||c 
(i.e. ^ab) is common to active and passive bands. In con- 
trast, the plateau in Cc/T vs. H at _ff^2 indicates that 
Hc2\\ab (i-G- £,c) is different between active and passive 
bands: i?,o||„h = , iJ' „„, = ^J" „ . The inter- 



2<a6Cc' c2\\ab 2iTC^bea' 

layer coherence length is expected to be shorter for the 
active 7-band with the dominant quasi- 2-dimensional dxy 
orbital character, thus limiting the main Hc2- In fact, the 
quantum oscillations indicate that the /3-band with the 
dyz and dzx orbital character has the largest interlayer 
dispersion.^ 

Next, we examine the polar angle 9 dependence 
of C4{9) under the condition of H/H^2{0) = 0.4 > 
Hth{^) / Hc2{0) ■ Figure 6 shows Cg/T at several polar 
angle 9, by rotating H conically as a function of </> and 
keeping 9 constant under the temperature T = 0.12 K. 



The field H at each 9 is adjusted so that the reduced 
magnetic field is maintained at H/Hc2{9) = 0.4, shown 
in Fig 5(a) as triangles. A non-sinusoidal 4-fold angular 
variation approximated as Ce{9, 0) = Co(0) + C4(^)/4((/)) 
is observed in the range of 80° < 9 < 90° . The oscillation 
amplitude C4{9) decreases steeply with decreasing 9 from 
90° and we cannot detect any angular (j) variation of Ce 
in the range of ^ < 80° . The results of the polar angle 6 
dependence of C4{9)/Cn and C4(6')/C4(90°) at T = 0.12 
K and H/Hc2 = 0.4 are summarized in Fig. 5(b). 

Now, we try to qualitatively estimate the evolution 
of Ci{9)/C4{^Q°) at zero temperature in a single band 
model. Since we keep the condition of exceeding the 
threshold field H,^ at 6 = 90°: > = 0.1, 

the expected result with nodes is qualitatively the same 
as that with gap minima. Thus we adopt the even- 
parity SC gap function A{k) = ^A^k^ky with verti- 
cal line nodes, instead of a realistic odd-parity A{k) = 
Z\o(sinafcj; -|- isinaky) for which there is no theoreti- 
cal results available so far. The calculation results of 
C4(6')/C4(90°), based on Refs. 26-28, with 2-dimensional 
(2D) cylindrical Fermi surface and 3-dimensional (3D) 
spherical Fermi surface are shown in Fig. 5(b). To in- 
terpret these calculation results, it is essential to remark 
that the Doppler shift 5oj = for k \\ H. In a 2D cylindri- 
cal Fermi surface case, since only the in-plane component 
of H gives a 4-fold angular variation for k \\ basal plane, 
C4(6')/C4(90°) is gradually suppressed, as H is tilted 
away from the basal plane and the in-plane component 
of H is decreased. This behavior agrees with the calcu- 
lation of point nodes in the basal plane with 3D spher- 
ical Fermi surface because the nodal k which involves 
the QP excitations by magnetic fields exists only in the 
basal plane. ^'''^^ In contrast, for vertical line nodes on 3D 
spherical Fermi surface, the behavior of C4(0)/C4(9O°) is 
more complicate and a 4-fold angular variation does not 
change by a small tilting H from the basal plane. The 
observed suppression of €4(9) shown in Fig. 5(b) is much 
more rapid than any of these single-band expectations. 

Therefore, we need to include the contribution of the 
SC gap structure in passive bands, in order to explain 
the polar angle 9 dependence of C4{9)/C4{90°). If the 
reduction of C4 is caused by compensation of anti-phase 
anisotropics between active and passive bands, the 9 de- 
pendence of C4 is expected to be rapidly suppressed when 
H < is satisfied, i.e., when both bands become rel- 
evant to anisotropic QP excitations. From Fig. 5(a), 
the applied field strength, described as the solid line of 
H/Hc2 = 0.4 where Ce(<?i) is measured, becomes smaller 
than -ffc2(^) below 9 = 81°. This characteristic angle cor- 
responds well with the observed threshold angle of steep 
suppression of C4{9) in Fig. 5(b). Thus we conclude that 
the main origin of steep suppression of C4 at low T and 
low H is the compensation by gap anisotropy with anti- 
phase between active and passive bands. Therefore the 
theoretical models (B)^^~^* is promising to explain the 
SC gap structure of Sr2Ru04; the SC gap in the passive 
bands most probably have the lines of gap minima or of 
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TABLE L The classification of typical gap structures in the 
superconducting state of Sr2Ru04 and consistency with ex- 
perimental results: Yes (Y) or No (N). 



#1 #2 #3 #4 #5 #e 



#7 



symmetry p p 

active band 
(7) 

passive band 





a. Ce(r) 


N 


N 


Y 


Y 


Y 


Y Y 


b. Cc{H) H\\ab 


N 


N 


N 


N 


N 


Y Y 


c. AC^h^T, 


N 


N 


N 


N 


N 


Y Y 


d. Ce(<^) 


N 


Y 


N 


Y 


N 


Y Y 


e. C^ifi) 


N 


N 


N 


N 


N 


N Y 


f. Scaling h||c 


N 


(Y) 


Y 


Y 


(Y) 


(Y) Y 


Summary 


N 


N 


N 


N 


N 


N Y 


Refs. 


14 


59,60 52-58 


65 


61-64 41-44 45-48 



zeros in the [110] direction. 



C. Superconducting gap structure of Sr2Ru04 

There have been many theoretical attempts to resolve 
the controversies on the SC mechanism of the spin-triplet 
superconductivity in Sr2Ru04 and on the nodal or node- 
like structure in its SC gap. The SC mechanisms me- 
diated by the anisotropy of the spin fluctuation^^"^^ or 
the orbital fluctuation,^^ due to the incommensurate AF 
fluctuation associated with a- and /3-bands, were pro- 
posed. Consequently, these models lead to the gap struc- 
ture which corresponds to phenomenologically proposed 
f^2_y2-wave state. '^^^^^ On the other hand, the spin- 
triplet superconductivity based on the third order per- 
turbation method for the repulsive Hubbard model for 
7-band, where the effective interaction beyond the spin 
fluctuation plays important roles, was discussed. The 
anisotropic p-wave state derived in this model has essen- 
tially the same gap structure as the one derived on the 
basis of the short-range ferromagnetic interaction. As 
other proposals, there are /-wave state with horizontal 
line nodes in the SC gap^^"^^ due to the interlayer pairing 
interaction and fxy-wave state with vertical line nodes. 

Furthermore, several theories, "'^"''^ taking the ODS 
into account, have been proposed. In these models, 
there are active and passive bands to the superconduc- 
tivity: the SC instability originates from the active band 
with a large gap amplitude; pair hopping across active to 
passive bands induces a small gap in the passive bands. 
The gap structure with horizontal lines of nodes^^"'*^ or 
strong in-plane anisotropy'*^"'** in the passive bands was 
proposed. 

Most of the proposed gap structure can be classified 
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FIG. 7: Schematic picture of the SC gap structure of Sr2Ru04 
on Fermi surfaces in fe-space. This figure is produced on the 
basis of the results in Ref. 45. 



in the five groups (#1— #5) with the single band models 
and the two groups (#6— #7) with the ODS models as 
summarized in Table I. In order to examine the consis- 
tency with the experiments, we discussed the items from 
(a) to (f). The discussions with respect to the items from 
(a) to (d) are included in our previous report.'^* Now, we 
try to examine them in the light of the present experi- 
mental results and select the appropriate gap structure to 
fully assign the gap structures in all bands of Sr2Ru04. 

(a) Cc/T in zero field shows T-linear dependence down 
to T/Tc = 1/30. This result requires the existence of 
lines of nodes, of zeros, or of deep gap minima in the 
SC gap in single band models, and in the case of ODS, 
the SC gap of the passive bands should have nodes, ze- 
ros, or gap minima. Thus #1 with an isotropic gap and 
#2 with shallow gap minima in single band models are 
inappropriate, (b) At low temperature, Cc/T for H \\ ab- 
plane reveals the typical field dependence with two kinds 
of the gap and supports the ODS models #6 and #7 
with the active band 7 and the passive bands a, (3. Next 
we address (c) the specific heat jump ACc/jnTc at Tc 
in zero field. The theoretical calculations of z1Co/7n7c 
with the weak coupling gives 1.43 for #1, 1.22 — 1.07 (for 
/imi„//imax - 1/2 - 1/4) for #2, and 0.95 for #3-#5 in 
single bands models, respectively, and thus overestimate 
the experimental result ziCo/7N?c = 0.73. On the other 
hand, the ODS models of #6 and #7, which have the gap 
structure of #2 for the active band 7, give good agree- 
ment: ACJ-f^Tc = (1.22 - 1.07) X 0.57 = 0.70 - 0.61, 
mainly originating from the active band 7 with the 170^ 
mimmMm (Z\inin/^max = 1/2—1/4). We note that a very 
deep gap minima in the active band will lead to smaller 
ACc/j-mTc than the observation. 

(d) For the field range 0.15 T < fioH < 1.2 T and 
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H 11 a6-plane, where the QPs in the active band 7 are 
the dominant source of in-plane anisotropy in Ce ((/)), we 
deduce the existence of a node or of a gap minimum along 
the [100] direction, because Ce takes a minimum with the 
4-fold oscillation. This result places severe constraints on 
the gap structure and requires the SC gap with the 4-fold 
anisotropy of #2 or #4 in active band, but not of #3. (e) 
Polar angle 6 dependence of the in-plane 4- fold oscillation 
C4^{9) gives an answer to the disappearance of the in- 
plane 4-fold oscillation of Ce((/') below 0.15 T. The issue 
was whether it originates solely from the gap minima 
of the active band 7 or mainly from the compensation 
for the field dependent QP excitations from the passive 
bands a,p. Steep suppression of C4{0) by tilting the 
polar angle 9 from the Ru02 plane is explained only by 
the compensation from the passive bands a,f3 with the 
ODS model #7. Finally, (f) the scaling of Ce{T,H) for 
H II c-axis confirms the existence of lines of nodes, of 
zeros, or of gap minima along the direction in the SC 
gap of both active and passive bands. 

To summarize consistencies with experiments on the 
items (a)— (f), we conclude that the ODS model #7 is 
suitable for description of the spin-triplet superconduc- 
tivity of Sr2Ru04. Figure 7 depicts the SC gap struc- 
ture of Sr2Ru04 on Fermi surfaces in fe-space, based 
on Ref. 45. If the essential order parameter is with the 
z{kx + iky) symmetry in the active band, the origins of 
the SC gap structure of #7 for each Fermi surface is ex- 
plained qualitatively as follows. The order parameter has 
nodes at (±7r, 0) and (0, ±7r) in the fc-space due to the 
periodicity and odd parity. For the active band 7, whose 
Fermi surface passes in the vicinity of these nodal points, 
the SC gap consequently has the lines of gap minima in 
the [±100] and [0±10] directions; thus the order parame- 
ter is described as d{k) =- zA^iswrnkx + isinafcj,).^'''^^ In 
the passive bands a and /3, the superconductivity with 
the same SC symmetry is coherently induced through the 
inter-band Coulomb interaction. However, the strong in- 
commensurate AF fluctuations by the nesting properties 
of the Fermi surfaces of a- and /3-bands, with the nest- 
ing vector in the diagonal directions, ^^"^^ give unfavor- 
able contributions to support this superconductivity with 
the z{kx + i^j/) symmetry. ^^""'^ Thus the SC gap in the 
passive bands have the lines of gap minima or zeros to 
the [±1±10] and [±1^10] directions to avoid the contri- 
butions of the incommensurate AF fluctuations, which 
would favor rf-wave. Therefore the experimental results 
are consistent with the models of #7, in which the order 



parameter with the z{k.j^ + \ky) symmetry in the active 
band originates from the pairing interaction which favors 
the odd parity orbital state, not from the ferromagnetic 
or antiferromagnetic spin fluctuations. 



IV. CONCLUSION 

In conclusion, we have obtained strong evidence for 
full determination of the superconducting gap structure 
of the spin-triplet superconductor Sr2Ru04 in all bands 
for the first time, based on the field-orientation depen- 
dence of specific heat at low temperature and the scaling 
of electronic specific heat in magnetic fields. We con- 
firmed the multi-band superconductivity with the active 
band 7 and the existence of a modulated superconduct- 
ing gap with the line of a minimum along the [100] di- 
rection in the active band. This gap structure in the ac- 
tive band is well described with p-wave order parameter 
d{k) = zAo{smakx + isinaky). Furthermore, in the pas- 
sive bands a and /3, the existence of the line of a minimum 
or of a zero in the induced superconducting gap to the 
[110] direction is experimentally indicated. This result 
suggests the incommensurate antiferromagnetic fluctua- 
tions at a- and /3-bands gives unfavorable contributions 
to this superconductivity. Among the superconducting 
mechanisms proposed so far for Sr2Ru04, the ones based 
on the vertex correction of Coulomb repulsion, other than 
spin fluctuations, best reproduce the experimental obser- 
vations. 
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